The linear growth rate of structure in Parametrized Post Friedmannian Universes 
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A possible solution to the dark energy problem is that Einstein's theory of general relativity 
is modified. A suite of models have been proposed that, in general, are unable to predict the 
correct amount of large scale structure in the distribution of galaxies or anisotropics in the Cosmic 
Microwave Background. It has been argued, however, that it should be possible to constrain a 
general class of theories of modified gravity by focusing on properties such as the growing mode, 
gravitational slip and the effective, time varying Newton's constant. We show that assuming certain 
physical requirements such as stability, metricity and gauge invariance, it is possible to come up with 
consistency conditions between these various parameters. In this paper we focus on theories which 
have, at most, 2"^^ derivatives in the metric variables and find restrictions that shed light on current 
and future experimental constraints without having to resort to a (as yet unknown) complete theory 
of modified gravity. We claim that future measurements of the growth of structure on small scales 
(i.e. from 1-200 h~^ Mpc) may lead to tight constraints on both dark energy and modified theories 
of gravity. 
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I. INTRODUCTION 

The dark energy problem, i.e. the possibility that 70% 
of the Universe seems to be permeated by an invisible 
fluid which behaves repulsively under gravity and does 
not cluster, has been the focus of research in cosmology 
for over decade. There are a host of proposals [l| and a 
battery of experiments are under way, or on the drawing 
board, to characterize the nature of this elusive source of 
energy 

In recent years, an alternative possibility has emerged, 
that Einstein's General theory of relativity is incorrect 
on cosmological scales and must be modified. Although 
the idea that General Relativity is incomplete has been 
around since the early 1960s p-lSl, there are now a num- 
ber of proposals for what this theory of modified grav- 
ity might be '10]. The Einstein-Hilbert action, Sg cx 
/ d^Xy/—gR (where g is the metric determinant and R 
is the scalar curvature of a metric gab) can be replaced 



by a more general form Sg 



-gF{R) where F 



is an appropriately chosen function of R illl, ll^; the 
dynamics of the gravitational field can emerge from a 
theory in higher dimensions such as one might encounter 
in brane worlds [l3| : a preferred reference frame may 
emerge from the spontaneous symmetry breaking of local 
Lorentz symmetry [l3-[l7|: the metric which satisfies the 
Einstein equation is not necessarily the one that defines 
geodesic motion [31 but is related to a second metric 
via additional fields (T9l - [2l| or connections j23 - [23 |: the 
Einstein-Hilbert action may be deformed by choosing as 
fundamental variables of gravity, SU{2) connections |25l - 



Many of these models have been successful in repro- 
ducing, for example, the observed relation between red- 
shift and luminosity distances from distant supernovae. 
They have, however, generally failed to reproduce the 
observed clustering of galaxies on large scales as well as 
the anisotropics in the Cosmic Microwave Background 
(CMB) unless the modified theory becomes effectively 
equivalent to general relativity (i.e. the Einstein-Hilbert 
action and a cosmological constant), e.g. [28l - [30| . The 
general problem that seems to plague most theories is 
an excess of power on the very largest scales which man- 
ifests itself through the Integrated Sachs Wolfe (ISW) 
effect and a mismatch between the normalization of the 
power spectrum of fluctuations on the largest and small- 
est scales. As yet, a truly compelling and viable model 
of modified theory of gravity has yet to be but forward 
which may resolve the dark energy problem. 

All is not lost, however, and progress can be made 
in learning about potential modifications to gravity by 
extracting phenomenogical properties that can be com- 
pared to observation- the " Parametrized Post Friedman- 
nian" approach j3l|. In this paper we focus on a key ob- 
servable characterizing the evolution of large scale struc- 
ture: the growing mode of gravitational collapse. 

The time evolution of the density field can be a sen- 
sitive probe of not only the expansion rate of the the 
Universe but also its matter content. In a fiat, matter 
dominated universe we have that 5m i the density con- 
trast of matter, evolves as 5m oc a where a is the scale 
factor of the Universe. We can parametrize deviations 
from this behaviour in terms of 7 [s^ - issj through 
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where Om is the fractional density of matter, % = ^ and 
overdots are derivatives with regards to conformal time, 
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T. For standard growth in the presence of a cosmological 
constant, one has that 7«6/lltoa very good approxi- 
mation. This is not true over a wide range of values for 
i^M- In fact, in Figure [T] we can see that 7 deviates from 
its early-universe asymptotic value as $7 a/ — ^ 0. A natu- 
ral question to ask is how 7 depends on different aspects 
of the Universe and how might use it to constrain dark 
energy and modifications to gravity. In this paper we will 
focus on a few of these properties. 
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FIG. 1: The solid line is the growth parameter, 7, for a ACDM 
universe, as a function of Qa/- For small values of f^A, 7 is well 
approximated by 6/11 (dashed line) but there are deviations 
as JIa grows; we find errors of 0.7%, 3.3% and 4.2% when 
= 0.7, 0.3 and 0.05. 

One important property of the Universe is the equation 
of state of dark energy, characterized by the constant (or 
function of time), w: 



Pe = wpE- 



(2) 



Pe and pE are the pressure and energy densities of the 
dark energy. The function w may be time varying and is 
related to the adiabatic speed of sound as 



w 



Cn =w 



m{l + w) 
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Another important property is gravitational slip, C, 
which is normally defined to be 



$ - * EE C$ 



(4) 



where we are taking a linearly perturbed metric in the 
conformal Newtonian gauge. 



ds^ = -a^(l + 2*)dr^ +a\l- 2$)dx 



(5) 



Such a parametrization has been advocated in a number 
of papers on modified gravity [sB-ES] and it has been 
shown that it can lead to a number of observational ef- 
fects. Albeit simple, and appealing, such a parametriza- 
tion of slip is not necessarily general and, as we shall see 
in the next section, necessarily implies other non-trivial 
modifications to the gravitational sector. Such modifi- 
cations are, in general, not explicitely acknowledged but 



may correspond to unexpected assumptions about any 
putative, underlying theory. Hence a more general as- 
sumption (at least within the context of 2"^^ order the- 
ories) would be that gravitational slip would depend on 
$ and $ (this is explained in more detail in section- Hl Dl 
and in [43]) 

Finally, we can define an effective Newton's constant in 
the relativistic Newton-Poisson equation 



X 



where Sx is the density contrast and 9x is the momentum 
of the cosmological fiuid X which has an equation of state 
wx- We can define the dimensionless function 
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eff 



(6) 



where G is the "bare" Newton constant. It then makes 
sense to try to constrain (7, w, ^ and p) in the hope that 
it may be possible to shed light on a possible theory of 
modified gravity. 

Although there alternative proposals [3, |4^ , a num- 
ber of groups have have pioneered the use of this simple 
parametrization of modified gravity (in terms of /^^ or 
both): in [s^, |33| it was argued that gravitational slip 
might be a generic prediction for modified theories of 
gravity, in (bsI . l39j it was shown that it would be pos- 
sible to constrain it through cross correlations of the 
CMB with galaxy surveys and in [i^ from the ISW ef- 
fect; weak lensing has been proposed as a possible route 
for constraining these parameters ,4'''" 49] with a tenta- 
tive detection of modification being proposed in [4l|, [s^] ■ 
Much is expected from applying these methods to future 
ambitious experiments that will map out the large scale 
structure of the Universe. Indeed constraints of General 
Relativity are a core element of the science that could be 
extracted from the Euclid experiment @ . 

Given that such an approach is phenomenological, the 
general attitude has been to leave these parameters com- 
pletely free. There is merit to such an approach in that 
one isn't restricting oneself to a particular theory and 
hence constraints will be general. It is true however that 
is is possible to idenitfy (reasonably general) consistency 
conditions for (7, w, C and /j,), contingent on specific 
physical assumptions. In this paper we state these as- 
sumptions and present restrictions on (7, wX and ji) . 
We shall use the formalism first proposed by one of us|43j 
which spells out how to build consistent modifications to 
gravity. 

This paper is structured as follows. In Section|TI]we re- 
cap the formalism presented in [4^ and relate it to the pa- 
rameters we wish to study phenomenologically. We dis- 
cuss how the consistency conditions reduce the freedom 
to choose arbitrary (7, p?). In Section |llT] we imple- 
ment the consistency conditions and find a relationship 
between the parameters by looking at the evolution equa- 
tion for the density contrast in matter for small wave- 
lengths. In doing so, we find find analytic expressions for 
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the relationships and briefly assess the range of scale to 
which they are applicable. In Section HVl we find analytic 
expressions for 7 to 2"'^ order for a general parametriza- 
tion which is consistent with the Parametrized Post New- 
tonian approximation on small scales. In Section|V]we dis- 
cuss the generality of the results and how they may be 
extended to other, more exotic models. 



II. THE FORMALISM 

We now summarize the formalism, the details of which 
can be found in the [i^: we present the field equations, 
the evolution equations for the fluid components and the 
consistency conditions for modifications to the field equa- 
tions. We shall further assume a spatially flat universe, 
but our results can be easily generalized to include cur- 
vature. 



A. The background cosmology 

As discussed in [4I], the background equations for 
any theory of gravity for which the metric is Friedman- 
Robertson- Walker (FRW) can be recast in the usual form 
used in GR. The Friedman equation simply reads 



(7) 



X 



In addition to the Friedman equation we also have the 
Raychaudhuri equation -2^+71^ = SnGa'^ J2x 



With the help of the Friedman equation, in a universe 
containing only pressureless matter and dark energy (as 
is approximately the case in the late universe) the Ray- 
chaudhuri equation may be rewritten as 
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The dark energy density pE and w, may be in gen- 
eral a function of additional degrees of freedom, the 
scale factor a or H. For example, for F{R) one gets 
PE = U^Fn - F)- ^Fr ~ ^Fr. But this ex- 
plicit dependence of pE (or of w) is irrelevant. One 
may always treat as a standard fluid with a time- 
varying equation of state w subject to energy conserva- 
tion PE + 37^(1 -I- w)pE = (but note that there may be 
additional field equations that determine the time depen- 
dence of w). In a universe containing only pressureless 
matter and dark energy, the energy conservation equa- 
tion for dark energy can be rewritten as 



(9) 



Our discussion above has one important consequence: 
that one cannot distinguish modifications of gravity from 
ordinary fluid dark energy using observables based on 
FRW alone. As discussed in 43], and further below, 
the situation changes drastically once we consider linear 
fluctuations. 



B. The field equations. 

The idea is to parameterize deviations from Einstein 
gravity at a linear level. Schematically we can write the 
modifled Einstein equations in the form 
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mod 



5Gau - 5Uau 



(10) 



Note that for a tensor F we use 6F to indicate a linear 
perturbation of F and we assume that SUat is made of 
the scalar metric perturbations and their derivatives. Let 
us also stress that the background tensor corresponding 
to Uab^ i-G. Uab vauishes. We assume that "normal mat- 
ter" (i.e. baryons, dark matter, neutrinos and photons) 
are contained in T^i, and that dark energy, or any non 
metric degrees of freedom that behave like dark energy 
(such as a scalar fleld- quintessence- or a dark fluid) , are 
contained in T^^. In this paper we will restrict ourselves 
to two fluids: T^i/ is the energy-momentum tensor for a 
pressureless fluid with density pM, density contrast dm 
and momentum 6]\i while is the energy momentum 
tensor of a fluid with density pde, density contrast Sde 
and Ode which can be characterized deflned in terms of 
(possibly time varying) equation of state and sound speed 
(we shall use the approach of [5l| to model a quintessence 
like fluid with a constant equation of state. 

The fleld equations can be rewritten in the following 
form 



-2fc2$ = SnGa^^pxlSx +Hl + wx)n9_ 

X 

+Ua + snug 

2($-t--H5') = 8TTGa^J2(^px +Px)ex + Ue 

X 

+ ^ A^Go'peTIe + Iup + Iw^U^ 
ar 6 3 



X 



(11) 
(12) 



$ - * 



Us (13) 



As advertised, the U terms contain modiflcations to grav- 
ity and we have used the notation from j43| : Ua = 
-a^U%, ViUe = -a'^U\, Up = a^U\ and [V^V^ - 
\V^5^j\Us = a2([/^ - \U\5^^). Further, we have pa- 
rameterized the dark energy pressure perturbation He = 
&Pe/pe as 



cl5E + 'i{cl~cl){l + w)'Hei 



(14) 



For adiabatic fluids, as is the case of radiation and CDM, 
Cs — Ca and we get He — c^Se- In general, however, 
Cs 7^ Ca and may in fact be a function of space as well as 
time. 



C. The fluid equations. 

It is convenient to define Ax = (5jc — 3(1 -|- wx)^ for 
X = M, DE. We then have that the equations of motion 
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for the fluids [Sll are: 



Am 



~n0M 

-9(l + w)-H(c2 



cD^E 



9e = {3cl-l)n9E 



{l + w)k^9E 

-cl)[<i> + n9E] 
1 



1 



w 



Ab + S* (15) 



D. The consistency conditions. 

In principle, one should be able to choose arbitrary 
combinations of metric perturbations to go in the tensor 
U. Yet in [43| it was argued that by assuming a set of 
general properties, it is possible to restrict the form of U. 
For the purpose of this paper, we shall choose the general 
theory of gravity to satisfy the following restrictions: 

1. The fundamental geometric degree of freedom is the 
metric. This encompasses most modified theories 
of gravity including the first order Palatini (torsion- 
less) formulations (where the connection F^j, is in- 
dependent of the metric at the level of the action) 
or purely affine theories [12] , provided a metric can 
be defined. 

2. The field equations are at most 2'^'^ order. This does 
restrict the class of acceptable theories (for exam- 
ple F{R) theories have generally higher derivatives) 
but it has become clear that it is these higher or- 
der terms that lead (again) to instabilities in the 
generation of large scale structure. 

3. The field equations are gauge-form invariant. 
Gauge form-invariance is the linearized version of 
the full diffeomorphism invariance of any gravita- 
tional theory with a manifold structure. It is the 
unbroken symmetry of the field equations under 
gauge transformations. After a gauge transforma- 
tion, the field equations retain their exact form : 
they are /orm-invariant (see [53l - l56| for further 
discussion). 

It is of course possible to relax some of these conditions 
and we will discuss how in the conclusions. 

Armed with these conditions we can construct U in 
Fourier space solely out of $ and $ such that: Ua — 
PA^, Ue = kB<i>, Up = fc^Ci* + fcCz* and f/s = 
Di^ + D2^/k. All operators above (e.g. A ) are di- 
mensionless. Terms such as $, 5* and ^E* (and higher 
derivatives) are forbidden while terms proportional to ^' 
are allowed but their coefficients vanish as a result of the 
Bianchi identities. This last fact is a consequence of the 
consistency conditions outlined above and the reader is 
referred to [11] for a more thorough treatment. 

We can rearrange (fTTj) in the form of ([6]) to immeadi- 
ately read off 

2 ^ , A + mkB 
M = 1 + ^ • 



where we let T-Lk = H/k. 

Furthermore we have that the shear field equation be- 
comes 



where we have let = Di and g = D2. 

We see that the gravitational slip $ — 5* will gener- 
ally depend on both <I> and its time derivative^ $, and 
not solely on Defining <E> — jppf'^ as has been in 
proposed in [s^ [s^ , we find that 



1ppf[^] = 



1 



-] _ /- 9 d In $ 

k (It 



g d In $ , , 



Thus unless g — 0, jppp is an explicit functional of 
introducing interesting enviromental dependance on the 
matter distribution. All parameterizations of the slip 
used so far, for which $ — \1/ oc have ignored this 
possibility which suggests that they were over simplistic 
(although see 1311') . 

As shown in [4J| , as a result of the Bianchi identities 
we have \/ aU\ = 0, leading to a series of restrictions on 
the coefficients: 



A = n 



2Hkg + 2k{2Ul + \)g-2HC, 



B 



A--g 

m 3^ 



2 



3 



Ci = -{B + 2HB)+2C 
k 



-A--[A + kB] 



Co 



-A 



(17) 



This means a consistent modification to the Einstein 
equations is uniquely determined by two arbitrary free 
functions, C(t, fc) and g{T,k). 

Finally we can combine the expression for A and B 
appropriately to find that g has a simple interpretation 
as a perturbation of the effective gravitational constant : 



2 . n 



(18) 



Hence the consistency conditions lead to an important 
relationship between a generalized form of the gravita- 
tional slip. In particular, if we consider time variations 
of the Newtonian potential, it is inconsistent to consider 
a restricted parametrization of the gravitational slip of 
the form <I> — = on all scales. [66|. 



E. Parameterizing and g 

We are not assuming any particular underlying the- 
ory of modified gravity and hence do not have a specific 
model for ^ and g. Our interest is in theories that may 
mimic the behaviour of dark energy so we expect devi- 
ations from Einstein gravity to emerge as VIe begins to 
diverge from 0. A simple assumption is to assume that 



5 



the gravitational slip is analytic at Qe = and Taylor 
expand it: 



We can do the same for g so that 



(19) 



(20) 



In Section UTT] we find how the growth depends on such 
a parametrization and, in particular, determine analytic 
expressions for 7. 

This way of parameterizing ^ and g has three major 
advantages : 

• It is in the spirit of the Parametrized Post- 
Newtonian (PPN) formalism where the PPN pa- 
rameters are isolated from the potentials which are 
dependent on the density profiles and thus the so- 
lutions; the role of the "potential" in this case is 
taken by ^IeIt) which depends on the background 
cosmology. 

• Expanding in powers offlE isolates the background 
effects of the dark energy from the genuine effects 
of the perturbations. In particular, the dark en- 
ergy relative density ^Iqe or the dark matter rela- 
tive density r^om would have no effect on the pa- 
rameters and gi. 

• This expansion makes mathematical sense for any 
analytic function, as the function fi£;(T) is always 
bounded to be < flE < 1, i.e. it is a naturaly 
small parameter. 

Note that we have dropped any k dependence from this 
parametrization. There are two ways that fc-dependence 
can enter, either relative to a fixed scale kg (which may be 
part of some theory of gravity) or relative to the Hubble 
scale H. If we wish to see how our results are affected 
by a scale dependence relative to the temporal changes 
introduced by the FRW background we can extend the 
parametrization to 



where, as above, we have 



^(1) 



(21) 



(22) 



and likewise with g'^^^ and g^^\ 

It turns out that, if we attempt to, on one hand gen- 
eralize our parametrization of ^ and g, but, on the other 
hand pin it down so as to be consistent with PPN method 
used on much smaller scales, we need to change our previ- 
ous approach. It is entirely possible that there are other 
scales in the system. Furthermore, as we show in in Ap- 
pendix 13 to leading order we may have g oc "H™ where 
m is negative. These can complicate the simple model 



we considered above. For example, consider the func- 
tion / — e^^ where ^ is a fixed scale. How does one 
expand this on small scales? One would be tempted to 
write / = e"^'- ^ 1 + iklik = 1 + iV. but this clearly 
makes no sense. The scale k was artificially introduced 
and leads to erroneous conclusions. 

In general we have a function C(r, k). Since C is dimen- 
sionless while r and k are not, the functional dependence 
on T and k must come in dimensionless combinations. 
It is convinient to exchange r with either r('H) or with 
t(S7£;). Thus the most general function ( will have the 
form ( — ({ilEjHkjiH, fc/fco) for constants £ and fco, and 
there may be additional dimensionless parameters enter- 
ing. We can thus isolate the leading-order dependence of 

on "Hfe and write 

C = a(r!£,fc)HI! (23) 

for a constant n. We expand g in a similar way as 

g^gLinE,k)nT (24) 

for a constant m. In Appendix [Xl we show that a consis- 
tent PPN limit fixes n — and to = — 1. 

Thus we arrive at our general expansion of C and g in 
the small scale limit, which is consistent with PPN: 



c 

9 

,(-1) 



9l^E + 92^E 



■ 93^% 



(25) 
(26) 



The parameters Q and gi may in principle be k- 
dependent, e.g. (1 oc (fc/fco)^ for a fixed scale ko and 
power index N. We shall not investigate this further 
in this work but we stress it as a possibility and note 
that our results would include these cases. In Section ITVl 
we conclude by presenting and analysing the resulting 
growth rate due to such a parametrization. 



III. THE GROWTH RATE ON SMALL SCALES 
FOR A SIMPLIFIED MODEL OF C AND g 

The definition of 7 originally arose when characteriz- 
ing the evolution of small scale density perturbations. 
We expect it to be particularly useful when characteriz- 
ing the growth of structure on small scales (by which we 
mean roughly between 1 and 200h~^ Mpc) as would be 
probed by galaxy redshift surveys (through redshift mea- 
surements of the pow er sp ectrum, for example, or redshift 
space distortions |57H59l |) and weak lensing surveys. 

In this section we focus on the behaviour of this sys- 
tem in the limit in which Hk = H/k <C 1, i.e. on 
scales deep inside the horizon. We can then assume that 
AoE — &DE — 0. This is true for ~ 0(1) or larger. 
Since in this paper we are concerned with modifications 
of gravity rather than the speed of sound we will leave 
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the full treatment of small for a future investigation. 
We shall, however, show the effect of small numerically 
futher below. 

In what follows we will present a modified evolution 
equation, find analytic expressions for the growth factor 
and compare to numerical results for the full system. 



A. Evolution of density perturbations 

Combining the fluid equations (jisp in one 2"^^ order 
equation, we find that Am obeys 



Km + nU Km - t^^mW V/^M = (27) 



with the damping coefficient modified by 
^[C,5] = 1 + 7771 TTTTif 



2(1-Hfe5) 

1 - Ukg + mlnM/2 

and the response term modified by. 

1 - C - 9B/2 



{l~'Hk9){l + mlnM/2-'Hk9) 



(28) 



(29) 



Specifying C and g completely fixes U and V. 

If we further take the provisional small scale limit 
^fc 1 (i-e. without assuming anything about C, and 
9) we find that 



nlr- 1 . , ^^MUk f 3Hk[l-C-9B/2] \ 

"t^'^^ ^ ' + ^0^^^^ + T^9 / 



(30) 



and 



V[C,.9] = 



{l-UugY 



The full small scale limit, including C and 9 is presented 
in appendix \K\ 



B. Analytic expressions for the growing mode 

We expect modifications of gravity to kick in when the 
expansion rate starts to deviate from matter domination. 
In this section we will work with the parametrization of 
C and g proposed in equations and We can 

immediately see from equations (|30|) and (|3T|) that the 
effects from 9 will only come in at order "Hfc . We shall also 
restrict ourselves to constant w and leave varying w for 
section ITVl In this section we shall present the derivation 
and result to first order vci^Ie and then present the result 
to second order vhVLe- 

The starting point is 



Am +'HA 



M 



M 



(l-Clf^B)A 



M 







Changing variables to In a and defining Aj\/ 
can rewrite this equation as 



oY , we 



Y" 



5 — iwVtE 



w - 



-CinM]Y = o 



where we have used the Raychauduri equation ([8]) and 
where we set ()' = -rr^- 

Changing variables to rig and using the O*'^ order fluid 
conservation equation, rewritten as fl'^ = —Sw^Im^e we 
find 



222 d'^Y 
3u; Qm^e ^0^2 



E 



w dY 

-rtM^^E [M2 - sfiE) - 5] -— 

Z as IE 

ll^Eii -w) + Cl^^M^^B]r = o (32) 



Equation ((32)) has a three regular singular points (at 
Ob = and Qe — ^) and can therefore be transformed 
into the hypergeometric equation. We wish to find its 
behaviour around CIe — and do so by expanding Y, 
Y = 1 + YiilE to find (to lowest order in He) 



1 



Ci 



w{5 — 6w) 
and hence 

1 - w + Ci 



A 



M 



1 



w{5 — 6w) 
Therefore to ©(I^b) 

l-w + Ci 



In A 



M 



In a 



w{5 — 6w) 



which we can use to find the logarithmic derivative of the 
growth factor / = din Am /dlna: 



(31) /=1 



3(1 



Ci) 



6w 



As stated above, we are parametrizing the growth factor 
using / 



fllj and so we have that 



7 = 70 = 



3( 1-w + Ci) 
5 — 6w 



(33) 



where the subsciript is due to the fact that this is the 
lowest order approximation. 

We can easily see that, for w = —1 and Ci = we 
retrieve 7 — 6/11 ~ 0.54545..., the approximation first 
proposed in [32| and subsequently rederived and advo- 
cated in [g^l and [s^. If we assume a more general (but 
constant) equation of state, we improve on the approxi- 
mation advocated in [33} . 

It is possible to further improve the approximation by 
going to next order in Qe- For this we need Y and ( to 
0{nl), i.e. Y = l + YinE + Y2nl and C = Cif^£ + C2f^|. 
. We now have that 



7 = 70 + Ji^E 



(34) 



7 



where 70 is as derived above and given in (p3)) and 
3w 



71 



[-3Yi - (2 - 3w)Y^ + 4^2] 



(35) 



where we have defined 

(1 - w){Ww'^ - 4w - 



l) + Ci(9w2 + 2w-2) 



2w2(l2w-5)(5-6w) 
Cf + w{5 - 6w)C2 



(36) 



2u;2(l2w-5)(5-6u;) 

In the case oi w = —1 (i.e. the cosmological constant) 
and C = we find that expression (I55|) reduces to 

15 
~ 2057 

As stated above assuming that g can be parametrized 
in the same way as independently of Hfc, we find that 
it does not aff'ect the growth of structure on small scales. 
The situation is of course different if we consider expand- 
ing 1 — /i^ = gHk in powers oOIe with coefficients which 
are Hk independent. 
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FIG. 2: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of Qm. The dashed curves are the 
numerical results for w — —1, —0.8, —0.6 and —0.4 in ascend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 



C. Comparison with numerical results 

We can solve P^ . (ITSI) and ([T7| to assess the qual- 
ity of this analytic approximation. We first restrict our- 
selves to a dark energy like fiuid with a large sound speed, 
cl ^ 0(1), (such as a quintessence model or most other 
field like models) and assume no modifications to grav- 
ity. We see a number of features in Figure [21 First of 
all, 7 is very clearly not independent of as has been 
generally assumed. In fact as VIm — > 0, 7 deviates sub- 
stantially from its asymptotic value at iljv/ = 1. Never- 
theless, we find that ((55)) is a good aproximation to the 
true behaviour. In Figure [2] we plot the true and approx- 
imate behaviours of 7 for w = —1, —0.8, —0.6 and —0.4: 
we find deviations of at most a ten percent for w — —0.4 
at =0.1. 

We may relax the condition — 1 substantially before 
the dark energy perturbations affect the growing mode 
in the density field. This is clearly illustrated in Figure [3] 
where, for w — —0.6, three different values of the sound 



speed are chosen: 
and = 0. For . 



= 5 
= 5 X 



X 10-4, 
10-4, 



ci = 10-^ cf = 10-'' 
the growth is still in- 



distinguishable from cj~l and only once the jean scale 
for the dark energy falls susbtantially below the cosmo- 
logical horizon is the effect noticeable. In particular our 
results hold for values Cg such that Cg > 10 ^ Hk, which 
approximately translates to >~ 3 x 10-*. 

We note in passing that the effects of the speed of 
sound without modifications of gravity, have been studied 
in (6l]-i63i] . In particular (63| have found fitting formulas 
for 7 which interpolate between Cg — and Cs ^ 0(1). 
However, those fitting formulas do not account for modi- 
fications of gravity and are in fact quite model dependent 
(they depend on the background cosmology). The cur- 
rent constraints on the speed of sound do not rule out 



small values [6J] and in fact they are consistent with 
Cs — 0. Since this work concerns the effects of modi- 
fications of gravity rather than effects coming from the 
speed of sound, however, we leave the case for small Cg 
for a future investigation. 
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FIG. 3: The growth parameter, 7, for a selection of dark 
energy models where w — —0.6 and the sound speed is chosen 



= 10' 



= 10"'' and c: = in 



descending order from the top of the figure, as a function of 

Let us now introduce modifications to gravity and as- 
sume a non-negligible gravitational slip. In Figures|l]and 
[5] we show how well the analytic approximation fares in 
comparison to different values of and C2 ( in these fig- 
ures we restrict ourselves to w = — 1 but the agreement 
between the numerical and approximate estimates of 7 
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is generic). In Figure |4] we restrict ourselves to a gravi- 
tational slip which is linear in fi^ and identify the two 
main effects. First of all, the quicker the onset of slip, the 
more effective the supression of growth due to the onset 
of dark energy- 70 increases with Ci. Furthermore, the 
dependence of 7 on fi^; , through 71 , changes sign so that 
for larger Ci, 7 becames smaller as f2j\/ decreases. 



0.75 
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0.65 
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0.55 
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"™ 

FIG. 4: The growth parameter, 7, for a selection of gravita- 
tional slip parameters of the form = (i^Ie, as a function of 
Qai. The dashed curves are the numerical results for ^1 = 0, 
0.2, 0.4 and 0.6 in ascending order and the corresponding an- 
alytic approximations are plotted in solid line. 

This last effect is further affected by C2- Indeed we 
find that the slope of 7 as function of Q m can greatly be 
affected by higher order terms in 



0.66 



0.64 



0.62 



0.6 




D. Intermediate and large scales 

While the 7 parametrization is particularly useful on 
small scales where terms dependent on Hk can be dis- 
carded, this isn't true once we look at horizon crossing, 
i.e. 7ik — 1- This regime will be of particular impor- 
tance for measurements of the cosmic microwave back- 
ground and using the Integrated Sachs- Wolfe effect to 
look for the presence of dark energy or modified gravity 
[ssl [sol . |46| . [65| . Let us now consider the parametrization 
presented in equations ([2T|) 

We expand eq. (|27p to include the first order term in 
Hk which gives 



n 



1 



ZHk^M (0) 



+nk 



(0) 



2-C 



(0) 



C 



(1) 







(37) 



where the functions C,'^'^\ ^'^''^nd g'-'^' are further ex- 
panded in powers of using equation (I22p . We apply 
the same techniques as in section llll Bl above. and expand 

7 as 

7-7("^+7(')Hfe 

where the coefficients 7*^°^ and 7^^^ are further expanded 



in powers of Q.e, i-e. 7 



(0) 



700 + 0{^e)- Carrying 



through the expansion we find that 



v(o) 



7(1) = _ 



3^^ + 0(17,) 
5 — bw 

27 3Co'i + 2Cii - 25, 

4 



01 



3w-l 



0{nE) 



(38) 
(39) 



There are now 3 free modified gravity constants: Coi, Cii 
and goi . Notice how the first order correction in Hk to 7 
only depends on the modified gravity parameters and is 
zero for wCDAI, i.e. corrections in Hk to 7 for wCDM 
come to second order. 

We find that scale dependent corrections in powers of 
Hk are always subdominant compared with corrections 
in powers of Qe, or corrections with respect to a fixed 
scale, e.g. ~ (ik)^ . They are thus effectively negligible. 
We find that at A: ~ O.OAhMpc^^ corrections in powers 
of Hk are around 1% at redshift z = 1 for = —0.6 
and become smaller at larger k (smaller scales), or lower 
redshift (^ ~ 0). Hence, it perfectly reasonable to discard 
scale dependent corrections which come in powers of Hk ■ 



0.58 
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FIG. 5: The growth parameter, 7, for a selection of gravi- 
tational slip parameters of the form = (^iQ,e + Caf^si as a 
function of Qm- The dashed curves are the numerical results 
for Ci = 0.2 and C2 = 0, 0.125, 0.2 and 0.4 and the corre- 
sponding analytic approximations are plotted in solid line. 



IV. GENERAL EVOLUTION OF ( AND g 

We now wish to address a more realistic expansion of 
the C and g proposed in the Section |TT] which is consistent 
with the PPN approximation on small scales, namely ([^5)) 
and (j26[) . Having in mind our findings in section llll Di on 
intermediate scales, we disregard any dependence of Am 
(and hence of 7) on Hk, and therefore to this order we 
can set Am = Sm- 
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We perform the calculation in steps. First, we solve the 
2nci order differential equation obeyed by 6 by applying a 
Taylor series expansion in f]^;, resulting in a set of coef- 
ficients {Yi, I2, Y3} which are functions of the expansion 
coefficients of w, U and V. Then we relate the a perturba- 
tive expansion coefficients for 6, namely {Yi, 12,^3}, to 
the 7 parameter. Finally we relate {11,1^2,^3} for gen- 
eral U and V to the specific case of our g(^CT)M model. 



A. Perturbative solution of the 5-equation 

We start with the equation for the matter density con- 
trast in the absence of DE perturbations, namely 



(40) 



Changing the independent variable from r to In a and the 
dependent variable from Sm to Y defined hy Y = 5m /o., 
we get 



Y" 



Y' 



y = 



(41) 



On small scales we may expand Y{k,T) = Y^^\k,T) + 
Y'-'^\k, T)'Hk + O('H^) (see appendixEl, where the func- 
tional coefficients F*^*'(/c,r) have no dependence on Hk 
but may still be /c or r dependent through combinations 
of the form £k or r/i (where £ is some scale, not neces- 
sarily the same scale for all such combinations). We are 
interested in the small-scale limit T-Lk — >■ 0, and we shall 
work with Y{k,T) = Y^^\£k,il,E) only. As discussed 
in section Ifff Dl above, at the scale of validity of the 7 
parameterization Hfe-corrections are always small and ir- 
relevant. Since time dependence only comes through t/£ 
for some scale i we may further exchange t/£ with a 
function of Qe, thus we let F(/c,r) = Y{k, Qe)- 

We now Taylor expand Y{k,nE) in powers of fi^;. To 
get 7 to 0(ri|;) we need to expand Y to 0{Qe) 



Y = 1 + Yini 



Y^ni 



Y.n 



(42) 



We then use the above expansion (|42|) into ([4T|) and 
match orders [g^]- To be able to do that we need to ex- 
pand the functions w, U and V. Since w always appears 
in the combination wVIe we only need it to 0{Q?e)- The 
functions U and V, however, are needed to 0(il|;). Thus 
we expand 



w 
U 
V 



Wo + WiflE + W2VI, 



-UiVlE 
-Vi^E 



V2^i. 



3"B 



(43) 
(44) 
(45) 



While Wi are constants, the Ui and Vi coefficients may be 
/c-dependent, for example Ui = Uoi{-^)^ for some index 



N and scale fco- Using the expansions ([42]), (|43]), p4)) 
and (|45p into (j4T|) and equating orders \x\VIe we find 



Yi 
Y2 



-1 + Wo + Vi - ^Ui 
Wo ( 6^0 - 5) 
1 



2wo(12wo — 5) [ 3 



Y. 



I - -U1 + V1+ 5wi 



4- Wo (-4 + 2Ui + 15wo - 18wi) 
1 f.. 2 



Yi 



3wo(18wo-5)l''^~3"^-''^+"^ 



V2 - -U2 - Vi - 4wi -I- 5w2 + 2wi {Ui - 6wi) 



+wo {-2Ui + 2U2 - 9wo + 42wi - 24w2) 



-I - -Ui+Vi- 9wo + lOwi 



Yi 



+2wo {2Ui + 27wo - 30wi) 



Yo 



(46) 



We notice that in general there are 9 initial coefficients 
appearing in (|40p that determine only 3 final coefficients 
Yi for the solution to (|40)) . 

Having found the coefficients Yi we proceed to relate 
them to 7. 



B. From 5 to 7 

We can now use the definition of the logarithmic 
change of the growth-rate / = and then get 7 from 
7 = In // In VIm ■ We find that 7 is expanded as 



7 = 70 + ll^E + 12^E 

where 



(47) 



70 
71 

72 



3woFi 



-70 



7o 



3 

2' 
11 



Yi 



1 



>7o 



:7o - 



L 2 



3wiyi + 6W0I2 
3 

:7i 



6 2" 2 

370 - 271] Yi 
-I- [—(1 — 6wo)7o + 6wi] Y2 + 9W0F3 



-- [6wi7o -I- 6w2 



(48) 



Given a set of coefficients {Yi, F2, Ya, wo, Wi, W2} we 
can get ^{^e)- Note that {¥1,12,^3} may be k- 
dependent, for example Yi = Yoi{-^)'^ for some index 
iV and scale fco- 

One important point is in order. What we have done 
so far is more general than the approach we discussed in 
the main part of the article. In particular the derivation 
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of the 7 coefEcients in this appendix would hold for any 
theory for which the density contrast obeys (PO)) . One 
such theory is DGP, even though strictly speaking DGP 
does not fit within our framework of the main part of the 
article. 

To connect the 7 coefficients above with our framework 
we must perform a third step : relate the lAi and V; 
coefficients with expansions of g and Q. 



The fit starts to diverge in the upturn of 7 for low Q.m 
a universal feature that becomes more pronounced for 
certain ranges of From inspecting these figures, it is 
clear that the behaviour of 7 near ^Ie — 1 — 51 a/ is telling 
something qualitative about C,. The parameter C,i drives 
the offset of 7 (i.e.70), C2 the slope (i.e. 71) and ^3 the 
curvature (i.e. 72). Measuring these three coefficents can 
give a direct handle on the time evolution of Q. 



C. Relating to the gCCDM model 

As shown in the appendix [Aj to be consistent with 
ultra-small scale quasistatic limit the functions C and g 
must have the form C, = + 0{T-Lk) and Hkg — g'"'^^ + 
0{Hk)- As discussed in the last part of section Hi E[ we 
further expand C'^^-' and g^^^'' as 



and 



(49) 



(50) 



respectively, where the coefficients may once again be k- 
dependent. To lowest order in Hk we find 



^1 = -^91 
3 



3 r 
2 [.93 



^3^7: 1.93 -32 + gl{gi 



1) + 25152] 



Vi = 251 -Ci 

V2 = 252 + (2 + 3wo)52-5iCi-C2 

V3 = 253 + 2(l + 3wo)5? + (3wi-3«;o-Ci)3? 

+5132(4 + 9u;o) - .92C1 - .91C2 - Cs (51) 



These expressions may then be used with (j46[) and 
to find the 7 coefficients. 



D. Comparison with Numerical Results 

The fitting formulas we have derived map theoretical 
properties of the gravitational field onto the observable, 
7. This allows us to circumvent the use of a full cosmo- 
logical perturbation code when trying to observationally 
constrain the ^ and 5 via the growth of structure. 
We have already seen how well such an approach fares 
for the simplified model we used in Section IlIII We now 
briefly show how the general framework fares- note that 
the model for 5 and C matches onto the PPN and we have 
found the expansion of 7 to second order in 51^;. Both of 
these properties should allow us to span a wide range of 
possible models. 

If we focus first on the gravitational slip, we can see 
in Figures HI [7] and [5] that the analytic fit works excep- 
tionally well, to within a few percent out to D,m =0.1. 
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FIG. 6: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of Qm. The dashed curves are the 
numerical results for ^1 — —0.6, —0.3, 0.3 and 0.6 in ascend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 
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FIG. 7: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of Qm. The dashed curves are the 
numerical results for ^2 ~ —0.6, —0.3, 0.3 and 0.6 in ascend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 

We find similar effect when looking at the time varying 
Newton constant, i.e. 5 (or ^^). Again there is a direct 
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FIG. 8: The growth parameter, 7, for a .selection of dark en- 
ergy models, as a function of Q,m ■ The dashed curves are the 
numerical results for fa = —0.6, —0.3, 0.3 and 0.6 in ascend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 



FIG. 10: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of flM ■ The dashed curves are the 
numerical results for 32 = —0.4, —0.2, 0.2 and 0.4 in descend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 



mapping between gi, g2 and 93 and 70, 71 and 72. The 
accuracy of the approximation breaks down for smaller 
values of JIm yet is still excellent in the range of interest 
for observational cosmology. For small values of D.e he 
accuracy is less than a percent and really only becomes 
large (of order 5 - 10%) for Om < 0.1. 
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FIG. 11: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of Q,m- The dashed curves are the 
numerical results for = —0.4, —0.2, 0.2 and 0.4 in descend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 



DISCUSSION 



FIG. 9: The growth parameter, 7, for a selection of dark en- 
ergy models, as a function of Q,m. The dashed curves are the 
numerical results for g\ = —0.4, —0.2, 0.2 and 0.4 in descend- 
ing order and the corresponding analytic approximations are 
plotted in solid line. 



Let us briefly recap what wc have done. The main 
point of this paper is that, when introducing modiflca- 
tions to gravity in linear perturbation theory, one must 
take into account the consistency conditions in the field 
equations. These necessarily lead to restrictions in the 
form of the modiflcations that can be introduced. Most 
notably, and within the context of second order theories, 
this means that if one wishes to include modifications to 
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the Newton-Poisson equation, then one cannot consider 
the simplified gravitational slip, 4* = and must in- 
clude an extra term such that = C$ + (g//c)(f> where 
^Ji = 1 — Geff/Go- If we wish to construct a proper 
"Parametrized Post Friedmanian" approach to modified 
gravity, any parameter we introduce must be independent 
of the environment or initial conditions in the perturba- 
tions. The only way to do this is to use the parametriza- 
tion we are advocating. To our knowledge, all attempts 
at studying cosmological deviations from general relativ- 
ity have ignored this and hence it is unclear what class of 
theories they map onto and which types of theories are 
being constrained. 

Having taken this point on board, we have found the 
expression for the growth parameter on small scales in 
terms of both the gravitational slip, C and the modified 
Newton constant, — 1 — ^H. Given a set of cosmo- 
logical constraints on 7 and its dependence on rij\/, it is 
now straightforward to calculate constraints on ^ and g. 
The growth parameter is given by equation (j47p which 
can be seen as a Taylor expansion in terms of 1 — i^M- 
The coefficients in this expansion, 70, 71 and 72 can be 
expressed in terms of the equation of state, w (see Equa- 
tion 143]) , C (see Equation [25]) and g (see Equation [26)) 
by using equations ([5T|) . followed by equations pS)) and 
finally equations pS)) . 

With these relations in hand, it is now possible to use 
cosmological observations to place constraints on the- 
ories of modified gravity. In this paper we have fo- 
cused on small scales (by which we mean between 1 and 
200h~^Mpc), scales that should be probed by redshift 
space distortion measurements, galaxy power spectra and 
weak lensing. Furthermore, we can now do this con- 
sistently, relating modifications in the growth rate with 
changes in the gravitational slip. This is of particular 
importance when considering weak lensing where obser- 
vations probe $ + It is also clear from our analysis 
that we have come up against the limitations of the 7 
parametrization: it is useful and effective on very small 
scales but not on scales comparable to the cosmological 
horizon. On those scales, one should be using the full set 
of field equations. We therefore do not advocate using 
our fitting formula to the growth on larger scales such as 
would be probed by the Integrated Sach- Wolfe. 

How general is this method? We have declared from 
the outset, the class of theories that we are considering. 
They must be metric, with 2"*^ order equations and sat- 
isfy gauge form-invariance. From what we have learnt 
about modifications of gravity, these seem a reasonable 
set of conditions to apply- they lead to theories which 
are less likely to be marred by gross instabilities either at 
the classical or quantum level. We should point out that 
all other attempts at developing such a parametrization 
have implicitely made these assumptions although have 
not necessarily done so self-consistently. It is possible to 
extend this analysis beyond the scope of these theories. 
If we are to go beyond 2"^^ order, one must include terms 
in $ or even higher. The Bianchi conditions will, again. 



impose a set of constraints on the coefficients of these 
terms and should allow a similar type of analysis. 

Two well studied theores are worth mentioning. F{R), 
F{R^'^ R^i,), etc, theories come with up-to four time 
derivatives in the field equations. Thus they do not fall 
directly within the methods of this paper but do under 
the general scheme outlined in [i^. In this case one 
would have to include terms involving <I>, $, ^ and ^ 
in to the Gqo and Goi Einstein equations, while the Gij 
equations would need <!> and ^E" in addition. Theories with 
higher derivatives are a subject that warrants further in- 
vestigation and have yet to be properly incorporated in 
any parametrized modifications of standard general rela- 
tivity. 

One other theory, studied extensively is the DGP the- 
ory In this case, only two time derivatives are 
present in the field equations and just like our frame, 
DGP contains two non-metric dynamical degrees of free- 
dom, which can be effectively written ad 6e and 9e- 
However, our framework cannot encompass DGP because 
DGP cannot be written as a generalized fluid as we have 
assumed of dark energy in this work (hence it is not a fail- 
ure of our use of SUab)- Nevertheless, our 7 parameteri- 
zation in powers of fi^ is still valid, and indeed needed. 
In the case of DGP we find that 



7 = — + -^^E - -^nl (52) 

16 5632 4096 ^ ^ ' 

gives an error on 7 around 5% at Um < 0.1, dropping 
to 2% at 0.2 and < 1% for larger values of flu- 

The error on the corresponding density contrast at those 
values of I^m is < 2%, ~ 1% and < 0.5% respectively. 
Notice how the coefficients are entirely fixed and do not 
depend on the only free parameter of the theory, namely 
the scale re- Rather Tc comes to play a role only through 

flE 
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Appendix A: The quasistatic limit : connecting to 
the Parametrized Post Newtonian (PPN) approach 

We reduce our equations to small scales and slow ex- 
pansions. To make contact with the PPN expansion, 
we will write the Einstein equations in powers of the 3- 
velocity v. The 3- velocity is related to the 4- velocity by 

= a~^w*. In this gauge, for scalar perturbations we get 
u' ^ a-'^W'e, so that = V'O. Hence 9 = -~k-'^W^v\ 
Letting v — \/ v^vi we get on dimensional grounds kd — v. 

The PPN order bookkeeping is k-^ = ()' ^ 0{v) and 
$ ^ ^ (5p ^ 0{v^). The same bookkeeping prescrip- 
tion holds in our case, and in addition, we also have 
Uk - 0{v) and Am Sm 0(0). 

We are now ready to find the small scale limit which 
is consistent with PPN. We start from the operators A, 
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B, Ci and Ca. Since = -jUHl + Siuf^fi) and n'j^ 
—3HkwV,M^E, and letting J — + wflE) we get 



2H,(.g'-C) + 2(2^2 + 



1 + |J 



S = 



1 + p 



Ci = -m 



1 



1 



l + P (l + fJ)2 
1 



2 + ml{l + 'AwVLe) 



1+|J^ ' (l + fj)2 

+27'H^l + wf^is + w2f^A/f^£;] \g' 



3H| 



^ 2C + 3-Hfe 



2.g" + 2-Hfe(2-3u'fl£;).g' 



and 

C2 = 3 
^ 6 



2Hk{9' -0 + 2{2Hl + \)g 



1 + |J 



where — > denotes taking only the lowest order terms that 
can contribute to the small scale limit. 

Now consider the Einstein equations. In the small scale 
limit we get 



2{^ + n^) = 8nGa^pe 

+n [-2C + 2g' + nk{i- 3wnE)g] $, 



(Al) 



(A2) 



6— + -H*) + 12-H($ + -H*) + 2k^{^ - *) = 
-3{EF + ER)^ + k^l^2(: + 3nk 2g" 
+HUi - GwUe + 9w^nE)g 

+2Hfe(2 - 3wf7£;)g'|* + 2 [-3HC + kg] (A3) 



and 



$ - * = + 

k 



(A4) 



As argued in Section we expand C in powers of Uk 
and write 

C = a(r!£,A:)HI! (A5) 

to leading order. We expand g in a similar way as 

g^gLinE,k)nT (A6) 

The goal now is to find the smallest powers m and n that 
can be consistent with the Einstein equations as Uk ^ 0. 

It is easily seen that (' = Ci^fe^^ for some function 
Ci (ils . k) which is found to be 



Cii = -77(1 + 3w^e)Cl - 'SwUm^e-^pt' 

Z UilE 



(A7) 



Similarly we have g' = and g" = 5l2'H^+^ and 

similar expressions to (jA7|) can be found for gLi and gL2- 
Consider again the Einstein equations and now keep 
on the lowest orders for each variable. For example $ is 
0(2) while {Ef + -Efl.)* = 0(4) and $ = 0(4) etc. We 
get 



-2A;^$ = 87rGaV5 + Gk'^Hk 



an 



n+l 



1 



(A8) 



+{1 - 3wnE)gLnk 



2CLni + 2gLin^+^ 

rn+l 



(A9) 



and 



5L2 + (2 - 3wf2_B)gLi 



~ iiwVlE + 9w''nE)gL 



(AlO) 



while the shear equation (jA4[) remains unchanged. 
Clearly the choice n = and m = —1 is consistent with 
all of the above equations. 
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Let us investigate whether smaUer numbers are pos- 
sible. Suppose that n < 0. Then the shear equation 



(|A4jl imphes (l'H^^ + ghHT^' = 0, hence if n < then 
m < — 1 and in particular m = n — 1. Note that this 
last relation includes the choice rt = as a special case. 
If on the other hand m < — 1 then the term gLhV^^' in 
the shear equation (jA4l) is of order less than two which 
forces automatically n < 0. Thus, without loss of gener- 
ality we may set m = n — 1. With this choice the Einstein 
equation (|A8p becomes 



Therefore if ji < we must have Ci^fc ~ \9l — and 
since both Cl and ql are independent of T-Lk, this forces 
Cl = gL — 0. Thus, the only consistent leading-order 
choice is n = and m = — 1. 

To summarize, a consistent small scale limit imposes 
the expansions 



C - CLinE,k) + oink) 

Tik 



(All) 
(A12) 



Note that there may be additional constraints on the 
fc-dependence of and gL- 

One source of worry is the "Hk^ term that persists in 
the shear equation (|A4[) . This is not a problem, how- 
ever. We may write the potentials in terms of the matter 
variables in a way that no ambiguity arises. We find 



Pm 



kHl-gL) 
l-a+5L (c 



M, 



(A13) 



gL 



i-gL 



2 1 - .gL 



'M, 



Pm 



1 



gL 



Om-H 



i-gL 



(A14) 
$, (A15) 



which are perfectly consistent equations. 

Finally, (jA14l) has a further interesting reduction (in 
this small scale limit). We replace 9m by —S/k'^, then use 
Sm — 'HfdM and finally use ()A13P to write $ — ^ppn'^ 
where 



1 



1 



ippN ^ — gL 



^ - Cl - gLil- gL)f 



+gL C.L - gL + 3wV,M^E 



dgL 
OQe 



(A16) 



while the measured Newton's constant on the Earth, Gn 
is 



G 



G 



N 



jppNi^-gL) 

Expanding (^e and gL in powers oi VIe we find 



7PPJV ~ I + (iflE 



and 



(A17) 



(A18) 



Gn 
~G 



hence 



1 + (.91 - <:i)^E 



(A19) 



N 



N 



-iw{gi - C,i)VLe'H 



(A20) 
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